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Abstract 

We show that the motion of a single hole in the infinite U Hubbard model with 
frustrated hopping leads to weak metallic antiferromagnetism of kinetic origin. 
An intimate relationship is demonstrated between the simplest versions of 
this problem in 1 and 2 dimensions, and two of the most subtle many body 
problems, namely the Heisenberg Bethe ring in 1-d and the 2-dimensional 
triangular lattice Heisenberg antiferromagnet. 
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The role of kinetic energy in the theory of magnetism is crucial: while virtual processes 
promote antiferromagnetism in insulators, as in the theory of superexchange, real ( or direct) 
kinetic processes usually promote ferromagnetism as a corollary of metallicity, as in the 
theory of double exchange [1,2]. The Nagaoka Thouless (NT) theorem [3,4] is of great 
importance in providing a rigorous mechanism for metallic ferromagnetism, notwithstanding 
the limitations of its context, that of a single hole in the limit of infinite repulsion. In this 
work we study the latter problem on certain 1-d and 2-d lattices with electronic frustration, 
a term defined by computing the sign of the hopping amplitudes around the smallest closed 
loop of a lattice, in complete parallel to the more familiar spin counterpart. If the sign is 
negative then the lattice is said to be electronically frustrated. The NT theorem applies 
only to the non frustrated cases, our focus is on the frustrated cases where not much is 
known reliably. We find surprisingly that in these cases the real kinetic processes promote 
antiferromagnetism in the metallic phase. 

Our immediate motivation for studying electronically frustrated lattices is to understand 
the physics of systems such as the recently found sodium cobalt oxide system Na x Co0 2 [5]. 
We present some analytical and numerical results for the one hole frustrated problem. Our 
main analytical tool is a novel reduction of the problem to an effective spin model, which 
yields insights into the physics. We have supplemented this with a numerical study using 
exact diagonalization of the effective spin problem, exploiting the symmetries of the clusters. 
Following the numerical work on the triangular lattice Heisenberg model (HM) by Bernu et 
al [6] we compute the exact eigenvalues in different total spin sectors for clusters of various 
sizes, and this leads to "towers of excitations" . These help one to distinguish between spin 
liquid states and various ordered states, something that variational studies cannot quite do 
[7]. Our findings are consistent with a three sublattice broken spin symmetric ground state, 
very similar to that of the triangular lattice HM [6,8] 

We study periodic clusters of the infinite U Hubbard model with L sites and N = L — 1 
particles with H = — Z)ij>*i,j c ^. (T c i| i a i where eg. a are Gutzwiller projected fermions with 
single occupancy constraint built into them. Due to the periodicity, we may define a wave 
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vector k for each wave function, and consider the action of if in a fixed k subspace. Let us 
locate the hole at site R in , and write a basis state la) = ct ct ...ct 10). A generic 
state may be written as ip(k,a) = -^J2r^ ik '^Tff\a), where T? is a (spatial ) translation 
operator 7>ct Tt = ct _ . Exploiting the translation invariance of the Hamiltonian, the 

— # 

matrix element of if in such a state is expressible in terms of an effective A; dependent 
operator 

mhfW) = ^tge^lTgc^c^Ja). (1) 

In this frame of reference, we hold the hole at a fixed site, as the entire set of spins flows 
past it. While this reduction is suggestive, it is not yet equivalent to finding an effective spin 
exchange representation, since the translation operators carry the complexity of fermionic 
negative signs. We can make considerable progress in specific problems as follows. 

1-dimension: Let us first consider the general 1-d case with arbitrary range of hopping 
t r . Let the hole reside at the Lth site, so that the L — 1 = N particles occupy sites 1 — > iV 
with some spin configuration \ai,..a N ) = c\ ai C2 a2 ..c' Ner \Q). Operating with terms in H 
shifting the hole to its left, we find tLj(-l) (L ^ ) cI iCT1 ..c]_ liCTj _ i c J t +1(Tj+l ..cj V)(Tjv c t Li(Tj |0). We 
next restore the hole to the Lth site from the jth site, so we need to translate by r = L — j 
units. This gives rise to negative signs that can be calculated readily, and the answer written 
down in terms of the permutation operator Pjj = | + 2Si.Sj, and spin translation operator 
T = Pi,2P2,3—Pn-i,n acting as a cyclic permutation on a lattice of length iV as 

H eff = J2(-l) r{L - r) (t r e-* kr P r , r+1 ...P 2 , 3 P 1 , 2 (Ty + he). 

r 

The sum is over the range of hopping, thus the spin problem is defined on a periodic ring 
of length iV rather than L — N + 1. In the simplest frustrated case we consider a model 
with nearest and second neighbor hopping t\ = t, t 2 = t', with both t,t' > 0. This (railroad 
trestle) lattice may be considered a strip of the triangular lattice. For this case, and with L 
odd, the effective Hamiltonian is 

H]f = te~ ik T + t'e- 2ik P lt2 {Tf + he, (2) 



the full spectrum is obtained by varying k. The hole has been eliminated, and we arrive at 
an impurity bond model residing on the deleted L — 1 site lattice. The first term of Eq(2) 
bodily translates all spin configurations, hence all spin states remain degenerate. The second 
term discriminates between the configurations through the single exchange term Pi j2 . 

The ground state of Eq(2) was found numerically for rings up to L = 19. L is taken to be 
odd to allow for the possibility of singlet ground states, and indeed in all cases, the ground 
state is a singlet [9]. The correlation functions are found by imposing pbc's with length 
N [10]. These alternate in sign and seem to decay as a power law for small t'/t, changing 
to a faster and possibly exponential decay at large t'/t. This behavior is reminiscent of 
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that of the HM with a second neighbor interaction, H h = Ji^2 n S n .S n+1 + J 2 J2 n S n -S n+ 2- 
Here it is known [11] that the Bethe point J 2 = has power law decay along with a 

logarithmic correction to the correlations so that asymptotically C(r) = J2i=i,NWi'^i+r) = 

i 

(— l) r { |7| + - B ( 1 °sj r D !? } + 0(l/r 2 ) [10]. The logarithmic term arises from the umklapp term, 
and B varies with J 2 within a power law phase persisting upto J 2 /Ji ~ .24 beyond which 
there is a gap in the spectrum. To compare with this behavior, we compute the structure 
function g(y,t') = Ef=i(-l) r C(r)- Plotted as a function of y = logiV in the power law 
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phase, it is expected to be a sum of two terms, one ~ y from A, and another ~ y^ from 
the logarithmic term B. For various values of t' we show g as function of y in Fig(l). We 
also display the corresponding structure functions on the same lattice sizes for the nn HM 
( J 2 = 0) and the Haldane Shastry (HS) model [12] which has B = 0. It is clear that the small 
t'/t cases are very similar to the Heisenberg model, whereas the case of t'/t large appear 
to be gapped. We note that in the limit t' — > + , the structure function g approaches the 
corresponding value for the HM. Remarkably enough, the individual correlation functions 
C(r) approach those for the nn HM to the available precision ( 6 decimal places)! Taking a 
closer look, we studied the wave functions for small clusters (upto L= 13). We found that the 
ground state of the impurity model, when translated as Y, r exp(iPr)(T) r | impurity gs > 
with an appropriate (Marshall sign) momentum P= (tt) for N/2 even (odd), becomes the 
exact ground state of the Heisenberg model. Based on these results we conjecture that the 
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t' — > + impurity model ground state after translation becomes the Bethe ground state of the 
Heisenberg ring. Finite t' seems related to further ranged exchange on a smaller scale than if 
( i.e. o(t')). We thus find that the impurity spin model may be regarded loosely as the HM, 
with a scaled down exchange J e ^ = t'/N to make meaningful comparisons of energetics. 

The Triangular Lattice In 2-dimensions the effective Hamiltonian Eq(l) can again be 
reduced to a spin operator. Let us imagine a two dimensional rhomboidal lattice with Ni 
columns and Mi rows ( L = NiMi) and the topology of a torus. We consider a particular 
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term where the hole has hopped to the site Ri — Si, and the translation is along Si so as 
to restore the hole to the fixed site Ri . It is enough to consider half the S^s corresponding 
to forward hops, the back hops contribute to the hermitean conjugate. The translation is 

— * 

expressible as a specific permutation of the L — l — (N\M\ — 1) variables Ri written in some 
specific order. This permutation can be decomposed into q cycles, each of length l m so that 
E m l m = L—l. The translation operator Tg. accumulates phase factors of (— \) lm ~ l from each 
cycle relative to the pure spin translation operator Tg , and one extra minus sign arises from 
the row containing the hole, so that the overall phase factor is -(-l)^™' !m_1 ' = (— l) L_Ci , 
and hence 

H k e ff = £{^(-1)^-^%- + h.c.}. (3) 

i 

For the triangular lattice we note that S = x gives q = Mi i.e. the number of rows, 
S = \x + ^-y gives q = Ni i.e. the number of columns. In the case of the third hop 
S = — 7}X + Ci depends upon the relative prime- ness of N\ and Mi but is trivial to 
compute for any given cluster once and for all. The translations T act only on the spin 
labels, and satisfy 7> ct Tj = ct for R ^ Ro + Si, and 7?cL^ TJ = cL ^ 
i.e. is a unit spin translator along the Si direction for all sites except the one nearest to the 
hole where it shifts by two units. 

To gain insight into the type of magnetic ordering in the triangular lattice we have per- 
formed exact diagonalization of small clusters [13]. We exploit the translation and rotation 
invariance by working with the deleted lattice in the subspace with S z = 0. We also use 
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time reversal invariance so that under a global transformation <jj — > —(Jj, its eigenvalues are 
k — il. We find that if (L — l)/2 is odd (even) the even spin states have k = —1(1), the 
odd spin states the opposite value. 

We contrast systems either supporting or frustrating 3-sublattice order. In particular, we 
choose systems of 9, 21 and 27 sites which support 3-sublattice order [13]. The 3 x 3-cluster 
displays high spatial symmetry which is reflected by a 7-fold degenerate ground state, one of 
which is a singlet which lies in the (k = 0, k — l)-sector. Due to its small size and geometry, 
in this cluster, every nearest neighbor (nn) is also a second neighbor and the third neighbor 
is the site itself. Thus, 3-sublattice order is forced by the boundary conditions in this 
cluster. We found a vanishing of the fluctuation of the operator Q = Y,<i,j>(Si 4 Sj — n i n j/4) 
in the singlet ground state. This curious result implies that within the tJ-model, both 
the contributions to the Hamiltonian share a common ground state for this cluster. In 
terms of the deleted lattice, the ground state of if e // is exactly the ground state of the nn 
HM. Next, we compute the spectra of the 21 and 27 site clusters. We found that the states 
corresponding to k = or Q* = ^-x were found to be lower in energy and nearly degenerate. 
This can be attributed to a unit cell tripling, as one expects in case of 3-sublattice order. 
We find in Fig(2) a systematic behavior of the excitations, in close parallel to those for the 
triangular lattice HM [6]. One can define a "moment of inertia" [6] / as the inverse of the 
slope of the line joining the bottoms of the different S tot towers of excitations. Our moment 
of inertia / ~ L 2 as compared to the Heisenberg value ~ L, again understandably in view 
of the extensivity of the latter model, and suggests J e // ~ t/L. A noteworthy feature is the 
striking "subgap" in each tower, this separates the ground state from the excitations that 
start out sparsely and then seem to form a continuum. A very similar feature in the HM 
has been identified with the magnons [6] with an energy scale uj = c\k\ ~ 1/yfL, and seems 
equally relevant here. In the 21-site cluster the states corresponding to three vectors ki, k 2 
and are all very close in energy. This can be understood in terms of a 1 st BZ diminished in 
size by 1/3 as these three momentum vectors are placed in equal distances from the corners 
of this new I s * BZ. In the 27-site cluster, the picture is very similar. Again, the k = or Q* 
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states are clearly separated from the rest of the spectrum. The family of first excited states 
is constituted by the momenta k±, k 3 and /c 4 , which are now vectors that lie exactly on the 
corners of the diminished 1 st BZ. In this cluster, a second family of excited states emerges, 
which is populated by states belonging to the k 2 subspace. This vector lies in the center of 
the neighboring diminished I s * BZ in y-direction of k . This point corresponds to a state 
of half the wave-number of the classical Neel state and is the least 'compatible' with the 
3-sublattice scenario, thus leading to the highest excitation energy. 

We also studied system of 15 and 21 sites, which frustrate the three sublattice order 
through their particular choice of boundary conditions [13]. We found that the energies of 
the frustrated clusters are considerably higher than for the unfrustrated ones. For 21 sites we 
can compare these directly: they are respectively —4.08577 and —4.18233, thereby showing 
the clear preference for three sublattice order. In plot (c) of Figure (2) we compare the tower 
of states obtained for this frustrated 21 site cluster. The regular features described above 
are now perturbed, the zero-momentum states no longer constitute the ground states of the 
different spin sectors and the quasi- degeneracy can no longer be identified. This scenario is 
similar to that obtained by a HM with an additional 4-site ring exchange term [14] . 

— # — # 

We also calculated the structure function S(q) = Y,i,j exp (iq ■ (f\ — r$)) < Si ■ Sj > / L. 
It is useful to identify v = \j S{Q*)/{L + 6) as the physical order parameter, with the same 
normalizations as in the HM [6]. The clusters compatible with 3-sublattice order contain the 

— * 

wave vector Q*-value, unlike the frustrated clusters. While the compatible clusters show a 
clear peak for S(Q*), the frustrated clusters still show a peak at the q- vector closest to Q*. 
Thus, 3-sublattice order is likely to be dominant in arbitrary systems and is not an artifact 
of a biased choice of boundary conditions. For the compatible systems we plot in the inset 
of Fig. (1) v vs. L~ x l 2 . and compare with results for the HM obtained in [6], normalizing 
to the maximum possible value of S(Q*). Extrapolation to infinite system size suggests a 
substantial finite intercept ( ~ 0.849). This shows strong similarities with the studies on the 
HM [6], and suggests Neel long range order originating solely from the motion of a single 
hole in a spin background. 
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In Fig (1 inset) we also plot the ground state energy. The small system size and high 
symmetry of the 9 site cluster leads to a fairly large ground state energy while the 21 and 27 
site clusters appear to show rapid convergence, we estimate E gs = —4.183 ± 0.005. On the 
square lattice, this swift convergence has also been observed by Poilblanc et al. [15]. Finally, 
the impurity model allows us to study the nature of the spin texture surrounding the hole 
[13]. In the (unfrustrated) hexagon surrounding the hole, we find substantial alternating 
antiferromagnetic order (Si ■ Sj) ~ —0.34 similar in magnitude and nature to that observed 
for a HM on a square lattice [16]. This local impurity-order rapidly transitions into three- 
sublattice antiferromagnetic correlations at greater distances from the hole. Thus, the hole 
can be seen as a moving impurity around which spins tend to line up antiferromagnetically. 

In conclusion, we have argued for kinetic antiferromagnetism on frustrated lattices. The 
preference for antiferromagnetism arises from the subtle phase dependence of the kinetic 
motion. Hole hopping mimics the effect of antiferromagnetic exchange energy, with the scale 
of a single impurity exchange bond per hole. The remarkable mapping to the Bethe ring in 
1-d and close similarity with the HM on triangular lattice highlight this subtle phenomenon. 
Our effect should be most prominent in situations where superexchange is negligible. At 
least qualitatively used, as in Ref [17], our findings can thus be interpreted as leading to 
weak antiferromagnetism with J e ff ~ +cx|i| in the case of the frustrated lattices, where x 
is the hole concentration. 
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FIG. 1. 1-d structure function g(y,t') of H 1 ' 1 ' for various i! compared to Heisenberg (H H ) and 
Haldane-Shastry (H HS ) model. In the limit t' — > 0, if' 1 *' correlations converge to those of . 
Inset: 2-d: Size dependence of the structure function on 2-d triangular lattice for Hubbard and 
Heisenberg [6] models denoted by solid triangles and open circles. The solid squares are the ground 
state energy of H e ff 
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momentum eigenvalue numbered in increasing magnitude of k, b: solid (open) symbols correspond 
to k = — 1 (k = 1) c: The frustrated 21 site cluster, zero momentum states (ko = 0) emphasized 
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Supplementary Plots 
3 sublattice-order preserving clusters 

Distinct colors correspond to different sublattices. 
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3 sublattice-order frustrating clusters 
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FIG. 6. The 21 site cluster frustrating 3-sublattice order. Check this by moving along the 
direction 1,2,3,... and placing 3 alternating classical Neel spins on consecutive sites, when arriving 
at the boundary a discontinuity arises. 
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FIG. 7. The 15 site frustrated cluster 
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